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Abstract 

The low-temperature series are calculated for the free energy, mag- 
netization and susceptibility in the Q-state Potts model on the square 
lattice, using the improved algorithm of the finite lattice method. The 
series are obtained to the order of z 41 for each of Q — 5 — 50, and the 
result of their Pade type analysis is compared with those of the large-Q 
expansion and the Monte Carlo simulations. 
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1 Introduction 



The finite-lattice method is a very efficient technique to generate the expan- 
sion series in statistical models ^] and in lattice gauge theory ||] . In this 
method the expansion series for the free energy density in the infinite- volume 
limit is given by an appropriate linear combination of the free energies on 
finite-size lattices. The partition function of the finite-size lattices, from 
which the free energy of the lattices is obtained, can be calculated using the 
transfer matrix formalism wm. This method avoids the problem involved in 
the graphical method, in which it is rather difficult to give the algorithm for 
listing all the diagrams that contribute to the relevant order of the series. 
The maximum order of the expansion series is determined by the maximum 
size of the lattices whose partition function can be calculated exactly to the 
order. 

In the case of the model whose spin-variable at each site takes many 
discrete values, however, one can only calculate the partition function of 
the finite-size lattices whose sizes are relatively small and it is difficult to 
obtain the series to such high order as in the Ising model, in which the spin 
variable takes only two values. Recently one of the authors || proposed 
an improved algorithm of the finite-lattice method, which enables one to 
obtain higher order terms than the standard finite-lattice method in the 
spin systems whose spin variables take a large or infinite number of discrete 
values. It was applied to the low-temperature expansion for the free energy 
of the solid-on-solid (SOS) model and it enabled him to extend the series 
to the order of u 23 § from the order of u 12 that was obtained using the 
standard finite-lattice method Q. 

In this paper we apply this improved finite lattice method to the calcu- 
lation of the low-temperature series for the free energy, magnetization and 
susceptibility in the Q-state Potts model on the square lattice. This system 
has the first order phase transition for Q > 5 |7|]. The values of the free 
energy, latent heat, and the gap of the specific heat and magnetization at 
the critical point are known exactly |||, while the exact values of the specific 
heat itself and the magnetic susceptibility at the critical point are not known 
yet. We calculate the series to order z 41 for each of Q = 5 — 50. We can 
extend the series from the ones calculated previously by Bhanot et al § and 
Briggs et al using the standard finite-lattice method, to substantially 
higher orders for large values of Q < 10, and our series for 11 < Q < 50 have 
the order that is about 3 times higher than the ones given by Kihara et al 
for general Q and zero field [12] or by Straley and Fisher [13] for general Q 
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and general field. We analyze the obtained series using the inhomogeneous 
differential approximants and their integration. It serves to check the va- 
lidity of the series obtained and of the method of the analysis used for the 
quantities whose values are known exactly. On the other hand, it makes a 
new prediction for the quantities whose exact values are not known. 

In section 2, we describe the improved algorithm of the finite-lattice 
method for the low-temperature expansion of the Q-state Potts model on 
the square lattice. The series obtained by the improved algorithm of the 
finite-lattice method are given in section 3. We present the analysis of the 
series in section 4. The result of our analysis is compared with those of 
the large-Q expansion and the Monte Carlo simulations. Our results are 
summarized in section 5. 

2 Algorithm 

Here we give the improved algorithm of the finite-lattice method for the 
low-temperature expansion of the Q-state Potts model on the square lat- 
tice, which enables us to obtain longer series than the standard finite-lattice 
method. Let us consider the two-dimensional L x x L y rectangular lattice 
Ao- The free energy density in the infinite- volume limit is given by 

/ = - lim -^-ln[Z(A ;Q)], (1) 

with 

Z(A ;Q)= ]T exp[-/3 £ (1 - S Si , aj ) - h^l - 6 Si , )] (2) 

{0<Si<Q-l} <i,j> i 

where the integer variable Si at each site i of Ao is restricted to < Si < 
Q — l. The low-temperature series is calculated with respect to the expansion 
parameter z = exp (— (3). We take the boundary condition that all the 
variables outside Ao are fixed to be zero. 

We consider the set {A} of all two-dimensional rectangular sub-lattices 
of Ao- The sub-lattice A is denoted by its size l x x l y and its position in Ao. 
We define a set of integers £ q that consist of and (q — 1) integers among 
the set of integers {1, 2, . . . , Q — 1}. We define the H of A and £ q as 

ff(A;£ 9 ) = -ln[Z(A;£,)]. (3) 
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In the calculation of the partition function Z(A;£ q ), the variable Sj at each 
site i inside A is restricted to be one of the elements of £ g , and all the variables 
outside A are fixed to be zero. We define W of A and £ g recursively as 

W(A-^ q )=H(A;Q- ]T W(A';Q). (4) 

a'ca, 

Here the prime in the summation implies that the W(A' = A; £ / = £ q ) 
should be excluded in taking the summation. We note that the H(A;£ q ) 
and W (A; £ q ) depend on the size of A but not on its position and depend on 
the value of q but not on the details of the elements of the £ q . So we denote 
them as H(l x x l y ; q) and W(l x x l y ; q), then we can rewrite Eq. (|j) as 

W(l x x ly-q) = H(l x x l y ;q) 



E 

l x <lx,l y <ly 

2<q<q 



(l x - 4 + l)(l y -l' y + l)[ q ,_\ ) W(l' x x l' y -q). (5) 



We know 



H(L X x L y ; Q) = {L x -l x + l){L y -l y + l)( Q q _\\ 

q<Q 

xW(l x x l y ;q). (6) 

Taking the infinite- volume limit we obtain 

/ = T lim ——H(L x xLy;Q) 

E ( Q ~l)w(l x xl y ;q). (7) 



lx,ly,q<Q 



The magnetization M and susceptibility \ are given from this free energy 
density, by 

and 

<9 2 f 
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In the practical calculation we introduce the variable x = exp (h) — 1 [14], 
then 

— I —I 

dh \h=o- dx \x=o ( > 

and 

d 2 I <9 2 5 I 

and it is enough to evaluate the expansion series for the free energy density 
to the second order in x in order to obtain the zero-field magnetization and 
susceptibility. 



In the standard cluster expansion of the free energy 15, 1£] for this 
model, a cluster is composed of polymers and each of the polymers consists 
of the excited sites that are connected by nearest-neighbor bonds. An integer 
value Si E {1, 2, . . . , Q — 1} is put to each excited site i of the polymer. For 
each cluster we can define the subset of the set {1, 2, . . . , Q — 1} so that 
each element of this subset is found to be put on at least one of the excited 
sites of all the polymers the cluster is composed of. Then we can assign to 
the cluster an integer q so that (q — 1) is the number of the elements of the 
subset. We can prove Q that the Taylor expansion of the W(l x xl y ;q) with 
respect to z includes the contribution from all the clusters of polymers in 
the standard cluster expansion that have q = q and that can be embedded 
into the l x x l y lattice but cannot be embedded into any of its rectangular 
sub-lattices. 

So the series expansion of the W(l x x l y ; q) starts from the order z n ^ x ' ly ' q ' 
with n(l x , l y , q) given in the following. It is enough to consider only for l x < 
l y . We first define k = [(q— l)/(l x l y )}, where we denote [p] as the maximum 
integer that is less than p. In the case of k = ( namely q — 1 < l x l y ), the 
cluster that contributes to the lowest-order term of the expansion series of 
the W(l x x l y ] q) consists of a single polymer. We know in this case that 



2(k + l y ) + q ~ 2 for q<l x + l y (12) 
l x + l y + 2q-2 for q > l x + l y . (13) 



n(l x ,l y ,q) 



Examples of the clusters that correspond to these two cases in Eqs. ( |12| ) 
and (|l3| ) are presented in Fig.l (a) and (b), respectively. In the figures the 
integers denote the excited spin variables at the sites and the solid lines 
represent the excited bonds. In the case of k > 1 ( namely q — 1 > l x l y ), 
the cluster that contributes to the lowest order term is composed of k + 1 
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polymers, among which a polymer consists of q' = q — 1 — kl x l y sites and 
each of the other k polymers consists of the l x x l y sites. Then if q' < l x 2 

k{l x (l y + 1) + (l x + l)l y } + 2q' + 2i + l 
n(l x ,l y ,q) = { for <»(*+!) ( 14 ) 

for </ > i(i + 1) (15) 

where i = [y^ ] . Examples of the clusters that correspond to these two 
cases in Eqs. (Ill) and (15) are presented in Fig. 2 (a) and (b), respectively. 



If q' > l x 2 

n(l x , l y , q) = k{l x (l y + 1) + (l x + l)ly} + 2q' + 2l x + j + 1 (16) 
where j = [(q' — l x )/l x ]- Example of the cluster that corresponds to Eq. (|i~6|) 



is presented in Fig. 2 (c). To obtain the expansion series to order z , we 
should take into account all the combinations of l x , l y and q that satisfy 
n(l x , l y , q) < N in the summation of Eq. ([?]) and should evaluate each of the 
W's to the order of z N . 

We should stress that the improved algorithm described above enables 
us to calculate the series for large Q that is substantially longer than the 
case when one would apply the standard finite-lattice method naively || [Tl]] . 
In the latter case one would define 

H(A;C Q ) = -ln[Z(A;C Q )} (17) 

and 

W(A; e ) = H(A; fo) - ^ W(A'; £ Q ) (18) 

A'cA 

for the rectangular l x x l y sub-lattice A. Here the variable s, takes all integer 
values < Si < Q — 1 for each site % € A in the calculation of the partition 
function Z(A; £q). Then the free energy density in the infinite-volume limit 
would be given by 

f=Y J W{l x xl y ;i Q ). (19) 

The series expansion of this W(l x x l y ; ^q) starts from the order of z 2 ( lx+l y\ 
The clusters that give the contribution to this order are the clusters consist- 
ing of a single polymer, an example of which is depicted in Fig. 3. So for the 
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series expansion to order z , one should take into account all the rectangu- 
lar lattices that satisfy 2(1 x + l y ) < N in the summation of Eq. (|l9|). The 
improved algorithm might appear to be more complicated than this stan- 
dard algorithm, since the former involves the summation with respect to q 
as well as the summation with respect to l x x l y . The CPU time and mem- 
ory for evaluating the partition functions of the finite-size lattices needed 
to obtain the series to the same order are, however, much smaller in the 
improved algorithm. In order to obtain the series to z 41 , for instance, we 
should calculate the partition functions up to the lattice-size of 10 x 10 either 
in the standard algorithm or in the improved algorithm. In the standard al- 
gorithm, the partition function with q = Q is needed for this maximum size 
of the lattice. On the other hand, in the improved algorithm the partition 
function with q = 2 is enough for this size of the lattice. 



3 Series 

Using the improved algorithm of the finite-lattice method described in the 
previous section we have calculated the low-temperature series for the free 
energy density /, magnetization M and susceptibility x m the- zero field 
for the Potts model on the square lattice to order z 41 for Q = 5 — 50. The 
obtained coefficients are listed in Table || and § only for Q = 20 and 50, where 
we give the series for z = exp (/) instead of the free energy density /, as well 
as the series for the magnetization and susceptibility. If the reader would 
like to know the coefficients for the other values of Q, he can get them from 
the authors by e-mail or from the authors' web-site. |jl0]| The calculation was 
performed by FACOM-VPP500 at KEK (Tsukuba) and FACOM-VP2600 at 
Kyoto University Data Processing Center. The maximum size of the used 
main memory was about 500Mbytes. We have checked that each of the 
W(l x x l y ;qys in eq. (^) starts from the correct order of z n<ylx,ly,q ^ described 
in the previous section. 

The series were obtained previously using the standard finite-lattice 
method by Bhanot et al || to order z 25 for Q = 8 and by Briggs et al 
ID to order z 39 , z 39 , z 39 , z 35 , z 31 and z 31 for Q = 5,6,7,8,9 and 10, re- 
spectively. The series were also calculated by Kihara et al\12\ to order z 16 



for general Q and zero field and by Straley and Fisher [|T3|] for general Q and 
general field to order z 13 using the graphical method. The series coefficients 
obtained here by us are consistent with the coefficients obtained by these 
people to their order. 
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We note that in the standard finite-lattice method the maximum order 
of the obtainable expansion series reduces when one increases the value of 
the Q, while the improved algorithm enables us to calculate the series to 
the same order for all the Q that are larger than or equal to some value. In 
fact, we can calculate the series to the order of u for an arbitrarily large 
Q(> 17) using the W(l x , l y \ g)'s that have been evaluated here to obtain the 
series to order u 41 for Q = 17, although we have stopped it at Q = 50. This 
can be understood by noticing from Eqs. (|T^)-(|T6|) that the series expansion 
of the W(l x , l y ] q) for arbitrary l x , l y and q with q > 17 starts from the order 
that is higher than z 41 and does not contribute to the expansion series to 
the order of z 41 . 



4 Analysis of the series 

The Q-state Potts model with Q > 4 has a first-order phase transition. 
For the first-order phase transition, some important quantities, such as the 
susceptibility and specific heat, have finite values at the critical point, in 
contrast with the fact that these values are infinite for the second-order 
phase transition. 

The infiniteness of the critical value for the second-order phase transition 
sometimes makes the analysis easier; if the physical quantity would have a 
strong singularity, we would be able to determine the critical exponent with 
a satisfactory accuracy by the differential approximations such as the Pade 
analysis. 

On the other hand, in the case of the first-order phase transition, we 
have to integrate the differential approximants to determine the critical am- 
plitude. This integration, however, has a subtle problem discussed below, 
and as a result, we must choose the approximants carefully. 

In analyzing the expansion series, we use homogeneous or inhomogeneous 



first order differential approximants [17], in which the approximants to a 
function f(z) satisfy 

Qi(z)f(z) + Q 2 (z)f(z) = R(z). (20) 

Here, Qi(z), Q 2 {z) and R(z) are polynomials of order Mi, M 2 , and M, 
respectively, which are determined so that Eq. (|20|) is satisfied to the order 
of z N with N the maximum order of the expansion series of the f(z). If R(z) 
is identically zero, the approximant is just the D-log Pade approximant. 
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To obtain the value of f(z) at the critical point z = z c , we must integrate 
the function f(z) = {—Q2(z)f(z) + R(z)}/Qi(z) from to z c . For the first 
order phase transition, the denominator Qi(z) of f'(z) has in most cases 
a zero close to the critical point. The zero point is usually slightly above 
the critical point on the real axis of z. It is sometimes, however, below the 
critical point, then the integral of f'(z) will diverge and we should abandon 
such an approximant. Even if the zero point is above the critical point, it is 
so close to the critical point that its subtle fluctuation makes the fluctuation 
of the critical value f(z c ) rather large. 

At first, in order to examine the accuracy of the analysis of the series 
described above, we apply it to the free energy density f c at the critical point 
(z c = 1/(1 + ^/Q)), the latent heat AU and the magnetization gap AM, 
whose exact values were obtained by Baxter ||. In Table || the estimates of 
the free energy density, latent heat and magnetization gap from the longest 
series (JV = 41) are presented for Q > 7, and their exact values are in 
parentheses. These estimates obtained by the analysis of the series may 
contain errors by two reasons; One of them is the statistical fluctuation 
among the approximants with the different orders M\ , Mi and M of the 
polynomials Qi(z) and Qi{z) in Eq. ( f20"| ) but with the same order N of the 
original series. About these statistical errors we find two common tendency; 
i) If the order M is too large, the convergence of the approximants is not 
so good, and if the orders M\ and M?, are too small, the integral of f'(z) 
diverges frequently, therefore we restrict the order M to be M = — 1 (D-log 
Pade approximant) or M = - 6 (inhomogeneous differential approximants) 
and the orders Mi and M2 to satisfy \M\ — M<i\ < 10. ii) For smaller 
Q, some approximants are finite at the critical point but apparently far 
from the exact value, then we must exclude them. We present only the 
statistical errors in Table ^. The other error is from the finiteness of the 
series. If there is discrepancy between an estimate for the finite series and 
the corresponding exact value beyond the statistical error for each N, we 
may say that the discrepancy is caused by the finiteness of the series and 
the error arising from the finiteness of the series has the order of magnitude 
of this discrepancy. 

Now, we present the result of the individual analysis for the free energy 
density, the latent heat and the magnetization gap. 

For the free energy density f c , both the inhomogeneous and homogeneous 
differential approximants give well-converged results, the latter converging 
more excellently. The data in Table |3| are from the latter. For Q > 30, 
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the estimates from the D-log Pade approximants agree with the exact value 
within one standard deviation. On the other hand, for Q < 20, the estimates 
are larger than the exact value, and do not agree with it within one standard 
deviation. As an example in Fig.4, we present a plot of the estimates for 
Q = 20 versus the number N of terms of the truncated series. The difference 
between the estimates and the exact value is about 0.01 percent while the 
fluctuation of the approximants is about 0.005 percent, the former is about 2 
times larger than the latter. The ratio of the difference between the estimate 
and the exact value to the statistical fluctuation becomes larger for smaller 
Q. Nevertheless, the figure also exhibits the tendency that the difference 
between the estimates and the exact value decreases with the increase of the 
number of terms. (We should mention that the errorbars for 35 < N < 38 
are relatively small. It is because, for these N, some approximants have 
divergent integral of f'(z) and the remaining fewer approximants happen to 
be very close.) Therefore, this difference can be expected to shrink for longer 
series. It would, however, be difficult to extrapolate the data to N — > oo to 
estimate the exact value, if it were not known, since we do not know what 
asymptotic behavior the approximants should follow with respect to N. 

For the latent heat, the inhomogeneous differential approximants give 
more convergent results than the D-log Pade approximants, in contrast with 
the case of the free energy density. Then, by using the inhomogeneous 
differential approximants, we can obtain about 60 estimates of the latent 
heat for each model with Q > 7. For Q > 40, the estimates coincide with 

the exact value within one standard deviation. For Q < 30, although the 
estimates do not agree with the exact value within one standard deviation, 
we find that their difference becomes smaller for larger number of terms of 
the expansion series and so we can say that the difference would disappear 
for long enough series, as in the previous case of the free energy density. 

For the magnetization gap, about 10 estimates obtained by the D-log 
Pade approximants converge well for each Q, while those obtained by the 
inhomogeneous differential approximants do not. Although it does not agree 
with the exact value within one standard deviation except for Q ~ 50, 
the estimates themselves are very close to the exact value; for example, in 
the case of Q = 20, the difference between the estimates and the exact 
value is about 0.015 percent. Further Fig. 5 shows that the differences would 
disappear for larger number of terms as well as in the case of the free energy 
density and the latent heat. 

Now using the same method as the above, we estimate the values of the 
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susceptibility and specific heat at the critical point, whose exact values are 
not known. 

The behavior of the differential approximants for the susceptibility is like 
that for the magnetization gap. The D-log Pade approximants for the sus- 
ceptibility show better convergence than the inhomogeneous approximants, 
although the errors of the estimates for the susceptibility are larger than 
those for the magnetization gap. The estimates from the longest series 
(N = 41) are listed in Table HI, and our estimates for 7 < Q < 10 are con- 



sistent with the estimates obtained by Briggs et a/[ll]. For Q = 10, 15, 
20, Monte Carlo estimates are given by Janke et al\23\. As their definition 
of the susceptibility is different from ours in the normalization factor, the 
results in Ref |23| should multiplied by a factor of {q/(q — l)} 2 to fit into our 
definition. The estimates obtained from Ref ||| are 3.777± 0.035 (Q = 10), 
0.7052 ± 0.0033 (Q = 15) and 0.30242 ± 0.00050 (Q = 20). Our corre- 
sponding estimates are 2.78 ± 0.14 (Q = 10), 0.657 ± 0.017 (Q = 15) and 
0.2938 ± 0.0036 (Q = 20). Although there are differences over one standard 
deviation, which are especially large for Q = 10, we may expect that the 
differences would disappear for larger Q. 

For the specific heat, the estimates obtained by the D-log Pade approxi- 
mants give well converging results, while the inhomogeneous differential ap- 
proximants do not. We list the estimates from the longest series (N = 41) 
by the D-log Pade approximants in Table || The errors of the estimates for 
the specific heat are smaller than those for the susceptibility. Nevertheless, 
our estimates for 7 < Q < 10 are inconsistent with the estimates obtained by 
Briggs et a/[|ll]] . Therefore we should check the validity of our approximants 
in some way. 

Fortunately, by the duality relation || , we can obtain the exact value of 
the gap of the specific heat as follows ||; 



In Table |, we also list the estimates of the gap of the specific heat, where the 
specific heat is calculated using the high-temperature expansion series 
obtained from the low-temperature series by the duality. It is clear that 
for Q > 11 there exists the gap of the specific heat and more over almost 
all the estimated values of the gap of the specific heat coincide with the 
exact values given by Eq. ( ^T]) within the standard deviation. Therefore, we 
convince the validity of our estimates for the specific heat. 
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Finally we compare our estimates of the specific heat at the critical point 
with the estimates obtained from the large-Q expansion series that was 
recently obtained by Bhattacharya et al |18}|, and with the ones obtained 



from Monte Carlo simulations |19], 2C, [2l|, [2^, 23]. In Table y, we show 



our estimates of the specific heat again and the estimates obtained from 
the large-Q expansion series and the Monte Carlo simulations. For Q = 
7, 8, 10, 15, 20, 30, the values of the estimates from the large-Q expansion 
are sited in Ref [18], the others are calculated by the regularized logarithmic 
Pade approximants given in Ref [18]. Among the Monte Carlo results those 
of Ref. ]23] have small statistical errors and they are consistent with the 
results of the large-Q expansion for Q = 10, 15, 20. Our estimates from the 
low-temperature series are systematically smaller than the estimates from 
the large-Q expansion. In Fig. 6 and 7, we present a plot of the estimates of 
the specific heat C~ for Q = 20 and Q = 50, respectively, versus the number 
N of the terms of the series truncated. We can see that the estimates are 
approaching the result of the large-Q expansion from below for each Q and 
that they approach faster for larger Q. Thus we can expect that, if we 
extrapolate the estimates in some way to N — > oo for each Q, we would 
obtain the estimates that may be more consistent with the result of the 
large-Q expansion, although we do not know how to make the extrapolation 
correctly, as mentioned above. 



5 Summary 

We obtained the low-temperature series for the free energy, magnetization 
and susceptibility of the Potts model on the square lattice. 

Using an improved algorithm of the finite lattice method, we extended 
the series from those given by the standard algorithm of the finite lattice 
method. Our improved algorithm is more efficient for the higher-state Potts 
models. 

Using the new series, we calculated the critical values of the free energy 
density, latent heat, magnetization gap, susceptibility and specific heat. 

The estimated values of the former three are very close to the known 
exact values. The values of the latter two, which are not known exactly, 
are obtained for various Q with high precision. Especially, the estimates of 
the specific heat are very close to the estimates from the large Q expansion 
series. 
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Table 1: The low-temperature expansion coefficients for the free energy, 
magnetization and susceptibility for the Q = 20 square lattice Potts model. 



n 


z 


M 


X 





1 


1 





1 











2 











3 











4 


19 


-20 


19 


5 











6 


38 


-80 


152 


7 


684 


-1440 


2736 


8 


-589 


1460 


-2280 


9 


4104 


-12960 


36936 


10 


34466 


-109520 


323912 


11 


-46512 


138240 


-199728 


12 


483398 


-2040420 


7678261 


13 


1614240 


-7086240 


30536496 


14 


-2010808 


4671200 


19059128 


15 


48638556 


-259282080 


1242232920 


16 


66750002 


-417704060 


2800285550 


17 


94851648 


-1155653280 


11160075600 


18 


4264583402 


-27837269920 


166147691080 


19 


2711283660 


—28967068800 


314891078616 


20 


36660260676 


-326575493500 


2700423490565 


21 


350530337736 


-2773117509120 


20580457476816 


22 


226902725452 


-3504413044080 


47276982916744 


23 


5653101156048 


-54127038428640 


481746412796832 


24 


28469510466547 


-273739383621540 


2536124771313858 


25 


42704598074544 


-626534197990560 


8388189200399472 


26 


700159265957050 


-7500005076811440 


75083308577493088 


27 


2423368661949468 


-28449325073826720 


325937254907403792 


28 


7731033540174246 


-110572911429835580 


1480698106281593115 


29 


78953345947533348 


-956758345892651520 


10938227296679669016 


30 


233171802367220678 


-3290695663943877920 


45003875885372649120 


31 


1207925854600942620 


-17857155839057961120 


248795624942616383064 


32 


8578891752427054593 


-118198071379878114580 


1550417616117914081584 


33 


26279052771743353332 


-427018427440815394080 


6633731834136735494256 


34 


169448419303276320478 


-2669246990084147723840 


397433094084062965431 12 


35 


934014456459287376408 


- 14635750107536528768640 


219466703191297256269272 


36 


3342530604057094050901 


-59864493412976212314860 


1016877625830611545531997 


37 


22227466452797269128528 


-378821839108831184988000 


6119224442354906298876816 


38 


105216422923458425776276 


- 1862780050724742025689040 


31557668972616720906032400 


39 


449173545929305443593316 


-8659232479851255144790560 


157936029967724328819066816 


40 


2809270008061904464069897 


-52184954811154800744150940 


921011351332383770376427924 


41 


12477877823064242776774212 


-246220110417995423090111520 


4633614244759932587057525712 
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Table 2: The low-temperature expansion coefficients for the free energy, 
magnetization and susceptibility for the Q = 50 square lattice Potts model. 



n 


2 


M 


X 





1 


1 





1 











2 











3 











4 


49 


-50 


49 


5 











6 


98 


-200 


392 


7 


4704 


-9600 


18816 


8 


-4459 


11150 


-20580 


9 


28224 


-86400 


254016 


10 


659246 


-2019800 


6003872 


11 


-1025472 


3225600 


-7018368 


12 


12281948 


-50077050 


193485271 


13 


89987520 


-370953600 


1510567296 


14 


-162228808 


567974000 


-886343752 


15 


3720812256 


— 19045411200 


92683335360 


16 


12381176822 


-67093148150 


376836950000 


17 


4469044608 


-89631139200 


1066326898560 


IS 


895089855002 


-5548225826800 


32926830187720 


19 


1892364280800 


-13708568736000 


107149526869056 


20 


13423422398316 


-111911089584250 


924768638863685 


21 


206844895493376 


-1536698755564800 


11106335245654656 


22 


478049132085952 


-4409226660382200 


42663913301510704 


23 


5667428458933248 


-51068413852233600 


452593396732200192 


24 


49346646734282527 


-435539546923202850 


3808046989807658388 


25 


172491504405475584 


-1825893233439542400 


19671214965561745152 


26 


1938967535168345710 


-19591446789689634600 


194639282941507768168 


27 


13089424692665334048 


- 135834659906560636800 


1410212431991986624512 


28 


63921872954358377616 


-750734284293452018450 


8836405252541575381155 


29 


628985731889425518048 


-7147786492756957612800 


80157906305935393997376 


30 


3839485035589391850158 


-45982543354085211648800 


551763520397297550084960 


31 


23463259463706802592160 


-302788074318283787164800 


3888163728504881482313664 


32 


203694326578125205781553 


-2595462149212992934427950 


32749085439406082880398944 


33 


1225877434229949054953952 


- 16606999340819591568259200 


224791274766846749241616896 


34 


8555010161977020587142598 


-120783311584253799401189600 


1691607730655033745533673792 


35 


66800901314888122476718848 


-948877081581250695129561600 


13374510959732980701033577152 


36 


417354805969242736779031351 


-6276444566058217793336109650 


93992501752971816579399237857 


37 


3081300147534452125687872768 


-47482109079734479233388368000 


725062658320307588405549126016 


38 


22520630415621757023630981256 


-353898756249751908192492010600 


5521526378626835675327811477000 


39 


148545447485043078418570354656 


-2445977768521719097911967382400 


39997172375744231349328962685056 


40 


1106109615352046378739786168307 


-18559451854582591017244861914850 


308483818208942470837183549866604 


41 


7848119276197068664163827052832 


-135219203614185545984625196732800 


2311827629935984338141853691434752 
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Table 3: Estimates of the free energy density f c , the latent heat AU and 
the magnetization gap AM . 



Q 




fc 


AU 


AM 


7 


-0. 


.0430048 ± 0.0000024 


0.35720 ± 0.00883 


0.7681 ± 0.0060 






(-0.0431119) 


(0.35328) 


(0.749565) 


8 


-0. 


.0392927 ± 0.0000055 


0.49058 ± 0.00193 


0.8080 ± 0.0056 






(-0.0393533) 


(0.48636) 


(0.799837) 


9 


-0. 


.0360255 ± 0.0000044 


0.60375 ± 0.00219 


0.8394 ± 0.0020 






(-0.0360805) 


(0.59967) 


(0.833261) 


10 


-0. 


0332047 ± 0.0000035 


0.69832 ± 0.00182 


0.8605 ± 0.0020 






(-0.0332365) 


(0.69605) 


(0.857107) 


11 


-0. 


.0307291 ± 0.0000024 


0.78047 ± 0.00057 


0.87746 ± 0.00072 






(-0.0307586) 


(0.77860) 


(0.875107) 


12 


-0. 


.0285722 ± 0.0000025 


0.85122 ± 0.00027 


0.89053 ± 0.00047 






(-0.0285900) 


(0.84994) 


(0.888878) 


13 


-0. 


.0266652 ± 0.0000023 


0.91316 ± 0.00029 


0.90114 ± 0.00032 






(-0.0266817) 


(0.91216) 


(0.899999) 


14 


-0. 


0249831 ± 0.0000021 


0.96756 ± 0.00012 


0.90989 ± 0.00025 






(-0.0249931) 


(0.96689) 


(0.909086) 


15 


-0. 


.0234812 ± 0.0000018 


1.015966 ± 0.000090 


0.92094 ± 0.00124 






(-0.0234908) 


(1.015414) 


(0.916664) 


16 


-0. 


.0221416 ± 0.0000018 


1.059128 ± 0.000067 


0.92351 ± 0.00018 






(-0.0221471) 


(1.058749) 


(0.923075) 


17 


-0. 


.0209338 ± 0.0000015 


1.098023 ± 0.000054 


0.92888 ± 0.000012 






(-0.0209393) 


(1.097700) 


(0.928570) 


18 


-0. 


0198439 ± 0.0000013 


1.133140 ±0.000041 


0.93358 ± 0.00012 






(-0.0198486) 


(1.132916) 


(0.933332) 


19 


-0. 


,0188560 ±0.0000011 


1.165129 ±0.000038 


0.937683 ± 0.000094 






(-0.0188594) 


(1.164923) 


(0.937499) 


20 


-0. 


,0179567 ± 0.0000010 


1.194298 ±0.000036 


0.941318 ± 0.000008 






(-0.0179586) 


(1.194155) 


(0.941176) 


30 


-0. 


0120422 ± 0.0000004 


1.390466 ±0.000011 


0.962987 ± 0.000006 






(-0.0120425) 


(1.390442) 


(0.962963) 


40 


-0. 


,0089727 ± 0.0000002 


1.498232 ± 0.000008 


0.972980 ± 0.000004 






(-0.0089727) 


(1.498224) 


(0.972930) 


50 


-0. 


0071096 ± 0.0000001 


1.567525 ± 0.000003 


0.978726 ± 0.000001 






(-0.0071096) 


(1.567522) 


(0.978723) 
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Table 4: Estimates of the susceptibility 



Q x{z c ) 



5 


154 ± 19 


6 


49.5 ± 4.1 


7 


19.0 ± 1.3 


8 


8.09 ± 0.29 


9 


4.55 ± 0.20 


10 


2.78 ± 0.14 


11 


1.900 ± 0.047 


12 


1.372 ± 0.026 


13 


1.036 ± 0.027 


14 


0.815 ± 0.019 


15 


0.657 ± 0.017 


16 


0.539 ± 0.016 


17 


0.4490 ± 0.0080 


18 


0.3844 ± 0.0061 


19 


0.3347 ± 0.0047 


20 


0.2938 ± 0.0036 


30 


0.1195 ± 0.0002 


40 


0.0697 ± 0.0002 


50 


0.0478 ± 0.00004 
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Table 5: Estimates of the specific heat. 



Q 


V 


V 


AC V 


(series) 


AC V 


(exact) 


7 


52.98 ± 0.61 


68.01 ± 0.47 


15.0 


± 


0.76 


0. 


.2234 


8 


28.47 ± 0.47 


30.17 ± 0.06 


1.70 


± 


0.48 


0. 


3099 


9 


20.03 ± 0.27 


20.87 ± 0.05 


0.84 


± 


0.27 


0. 


3841 


10 


15.579 ± 0.083 


16.294 ± 0.034 


0.715 


± 


0.090 


0. 


.4476 


11 


12.856 ± 0.061 


13.493 ± 0.034 


0.637 


± 


0.070 


0. 


.5021 


12 


10.985 ± 0.050 


11.622 ± 0.055 


0.636 


± 


0.074 


0. 


.5492 


13 


9.598 ± 0.032 


10.258 ± 0.062 


0.659 


± 


0.070 


0. 


5901 


14 


8.548 ± 0.027 


9.204 ± 0.020 


0.656 


± 


0.034 


0. 


6260 


15 


7.721 ± 0.028 


8.388 ± 0.004 


0.666 


± 


0.029 


0. 


.6576 


16 


7.044 ± 0.011 


7.739 ± 0.002 


0.695 


± 


0.011 


0. 


.6856 


17 


6.4908 ± 0.0086 


7.2061 ± 0.0041 


0.715 


± 


0.010 


0. 


7106 


18 


6.0251 ± 0.0074 


6.7637 ± 0.0095 


0.739 


± 


0.012 


0. 


7330 


19 


5.6317 ± 0.0061 


6.3859 ± 0.0098 


0.754 


± 


0.012 


0. 


.7532 


20 


5.2913 ± 0.0053 


6.0645 ± 0.0069 


0.773 


± 


0.009 


0. 


.7714 


30 


3.4012 ± 0.0010 


4.2919 ± 0.0051 


0.8906 


± 


0.0052 


0. 


.8861 


40 


2.5864 ± 0.0006 


3.5303 ± 0.0019 


0.9439 


± 


0.0020 


0. 


9393 


50 


2.1218 ± 0.0004 


3.0936 ± 0.0007 


0.9717 


± 


0.0008 


0. 


.9668 
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Table 6: Estimates of the specific heat: low-temperature series vs. large Q 
series and Monte Carlo simulations. 



Q low-temperature series 



large Q series 



Monte Carlo simulations 



52.98 ± 0.61 



28.47 ± 0.47 



69.6 ± 0.5 



36.9 ± 0.2 



9 


20.03 ± 0.27 


24.1 ± 0.2 


10 


15.579 ± 0.083 


17.98 ± 0.02 


11 


12.856 ± 0.061 


14.24 ± 0.05 


12 


10.985 ± 0.050 


11.84 ± 0.02 


13 


9.598 ± 0.032 


10.17 ± 0.01 


14 


8.548 ± 0.027 


8.939 ± 0.009 


15 


7.721 ± 0.028 


7.999 ± 0.003 


16 


7.044 ± 0.011 


7.248 ± 0.004 


17 


6.4908 ± 0.0086 


6.643 ± 0.003 


18 


6.0251 ± 0.0074 


6.142 ± 0.002 


19 


5.6317 ± 0.0061 


5.720 ± 0.001 


20 


5.2913 ± 0.0053 


5.3612 ± 0.0004 


30 


3.4012 ± 0.0010 


3.41294 ± 0.00005 


40 


2.5864 ± 0.0006 


2.58986 ± 0.00002 


50 


2.1218 ± 0.0004 


2.123198 ± 0.000006 



50 ± 10 p] 
47.5 ± 2.4 EI | 
44.4 ± 2.2 §2| 



12.7 ± 0.3 |19] 
17.81 ± 0.10 12I 



8.004 ± 0.019 §3| 



5.2 ± 0.2 @ 
5.3612 ± 0.0055 || 
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Figure Captions 



Fig. 1: Examples of the cluster consisting of a single polymer ( k = 0) 
that contributes to the lowest-order term of the expansion series for (a) 
W(l x = 4,l y =5;q = 4) and (b) W(l x = 4,l y = 5;q = 11). 

Fig. 2: Examples of the cluster consisting of two polymers ( k = 1) 
that contributes to the lowest-order term of the expansion series for (a) 
W(l x = 3,l y = 4;q = 18), (b) W(l x = 3,l y = 4;q = 21) and (c)W(l x = 
3,l y = 4;q = 24). 

Fig. 3: Examples of the cluster consisting of a single polymer that 
contributes to the lowest-order term of the expansion series for W(l x = 
±,l y = 5;£q). 

Fig. 4: The estimates of the free energy density f c for Q = 20 versus the 
number of terms of the series. The solid line shows the exact value. 

Fig. 5: The estimates of the magnetization gap AM for Q = 20 versus 
the number of terms of the series. The solid line shows the exact value. 

Fig. 6: The estimates of the specific heat C~ for Q = 20 versus the 
number of terms of the series. The solid line shows the estimates obtained 
from the large-Q expansion series by using the regularized logarithmic Pade 
approximants. 

Fig. 7: The estimates of the specific heat C~ for Q = 50 versus the 
number of terms of the series. The solid line shows the estimates obtained 
from the large-Q expansion series by using the regularized logarithmic Pade 
approximants. 
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Figure 1: Examples of the cluster consisting of a single polymer ( k = 0) 
that contributes to the lowest-order term of the expansion series for (a) 
W(l x =4,l y =5;q = A) and (b) W(l x = 4,l y = 5;q= 11). 
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Figure 2: Examples of the cluster consisting of two polymers ( k = 1) 
that contributes to the lowest-order term of the expansion series for (a) 
W(l x = 3,l y =A;q= 18), (b) W(l x = 3,l y = 4;q = 21) and (c)W{l x = 
3,l y =4;g = 24). 
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Figure 3: Examples of the cluster consisting of a single polymer that 
contributes to the lowest-order term of the expansion series for W(l x = 
4,l y = 5;£q)- 
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Figure 4: The estimates of the free energy density f c for Q = 20 versus the 
number of terms of the series. The solid line shows the exact value. 
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Figure 5: The estimates of the magnetization gap AM for Q = 20 versus 
the number of terms of the series. The solid line shows the exact value. 
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Figure 6: The estimates of the specific heat C~ for Q = 20 versus the 
number of terms of the series. The solid line shows the estimates obtained 
from the large Q expansion series by using the regularized logarithmic Pade 
approximants. 
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Figure 7: The estimates of the specific heat C~ for Q = 50 versus the 
number of terms of the series. The solid line shows the estimates obtained 
from the large Q expansion series by using the regularized logarithmic Pade 
approximants. 
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